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Abstract 



O ! We calculate the spin susceptibility at and below T c for a d x 2_ v 2-w&ve super- 
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Q^ I the impurity scattering and the Coulomb correlations act to maintain peaks 



conductor with resonant impurity scattering and Coulomb correlations. Both 



in the spin susceptibility, as a function of momentum, at the Brillouin zone 
edge. These peaks would otherwise be suppressed by the superconducting 
q , gap. The predicted amount of suppression of the spin susceptibility in the 
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superconducting state compared to the normal state is in qualitative agree- 
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rN | ment with results from recent magnetic neutron scattering experiments on 



Lai.86Sro.i4Cu04 for momentum values at the zone edge and along the zone 

diagonal. The predicted peak widths in the superconducting state, however, 

are narrower than those in the normal state, a narrowing which has not been 

observed experimentally. 
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I. INTRODUCTION 

Recent magnetic neutron scattering experiments on La 186 Sr .i4CuO4 by Mason and co- 
workerstl show an isotropic but incomplete suppression of the scattering intensity below T c 
for small energy transfers and momentum transfers near the (tc, tc) point of the Brillouin 
zone. In the normal state, the magnetic neutron scattering intensity, as a function of mo- 
mentum, exhibits four peaks displaced slightly along the edges of the Brillouin zone away 
from the (tc, tc) point at the corner of the zone, as shown in Fig. [H The experimental results 
show a suppression of these peaks, compared to their magnitude at T c , of about 60% at a 
temperature of 0.13T C . 

For an isotropic s-wave superconductor we would expect an isotropic suppression of this 
scattering intensity as the temperature drops below T c . But for a clean s-wave state at low 
temperatures the intensity should die away exponentially with temperature. The moderate 
suppression shown by experiment suggests that if s-wave superconductivity is present, it is 
gapless and isotropic. 

Alternatively, a superconductor with d x 2^ y 2 symmetry is always gapless, so an incomplete 
suppression of the scattering intensity is to be expected. However, for a clean ci-wave 
superconductor it is expected that this suppression will not be isotropic. LuH and Zha and 
co-workersa have pointed out that in the superconducting state the four normal state peaks 
should be suppressed, leaving only the response due to particle-hole excitations associated 
with the nodes. This response is displaced slightly away from (tc, tc) along the zone diagonals, 
giving a low temperature structure which appears to have undergone a 45 degree rotation 
about (tc,tc). The failure to observe this rotation has raised questions regarding whether a 
d x 2_ y 2 description of the superconducting state of La2_ a; Sr a ;Cu04 is viable. 

Now, for d x 2_ y 2 pairing it is known that resonant impurity scattering leads to important 
effects for both the penetration deptho and the NMR responseB In addition, the effects of 
the Coulomb enhanced spin fluctuations are believed to play an essential role in the cuprates. 
As previously discussed,™! when the effects of Coulomb correlations on the magnetic neutron 




FIG. 1. Map of reciprocal space probed by the neutron scattering experiments. Dotted lines 
represent the Brillouin zone edges near the (ir, it) point. The dashed line indicates the zone diag- 
onal. Filled circles show the locations of neutron scattering intensity peaks in the normal state. 
Open circles represent the expected positions for scattering intensity peaks for a clean d-wave 
superconductor well below T c . The solid line shows the direction of the momentum scans shown 
in Fig. H and Figs. |lO| - |T^ . 



scattering are included within the RPA approximation, the d x 2_ y 2 results can be brought 
into closer agreement with the experimental observations. 

Here we investigate both of these effects by calculating the neutron scattering for a model 
of a dirty d x 2_ y 2 -wave superconductor with strong Coulomb correlations. Specifically, we 
consider a Hubbard model on a two-dimensional square lattice with a near-neighbor hopping 
t and an on-site Coulomb interaction U. To this we add impurity scattering in the form of 
a dilute random array of zero-range scattering potentials V. 

The magnetic neutron scattering intensity is given by a structure factor >S(q, uo) = [n(u>) + 
l]x"(q, u>), which is proportional to the imaginary part of the spin susceptibility x(q 5 w). We 
calculate x(q, u>) in two stages. First, in Sec. ||, we examine the effects of the impurities 



on x(q, u), including self-energy as well as vertex corrections. Then, in Sec. |HI| , the strong 



spin susceptibility enhancement effects of the Coulomb interaction are included within an 



RPA-type approximation. Section [iy| contains our conclusions. 

Before beginning, it is useful to review the behavior expected for a system described by 
a two-dimensional tight-binding band which undergoes a transition from a normal state to a 
d x 2_ y 2 BCS superconducting state in the absence of impurity scattering and strong Coulomb 
correlations. As discussed in Littlewood et al.p the four normal state peaks shown in Fig. [I] 
arise due to a combination of the weak requirements for favorable nesting in two dimensions 
combined with umklapp scattering processes. As a d-wave gaps opens up below T c the 
thermally available scattering states are restricted to those near the nodes of the <i-wave 
gap function and these peaks are suppressed. The favorable low-energy momentum transfers 
are then restricted to those which represent node-to-node scattering, and the corresponding 
wavevectors lie along the zone diagonal directions. Thus at low temperatures and energies 
well below 2A(0), peaks in the neutron scattering should be observed to lie along the zone 
diagonals.™ 

Here we find that impurity scattering and Coulomb correlations both act to keep the 
scattering intensity peaks at their normal state positions along the Brillouin zone edge. 
For reasonable parameters, the amount of suppression of the scattering intensity in the 
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FIG. 2. Lowest order bubble diagram representing the spin susceptibility probed in magnetic 
neutron scattering experiments. 

dirty-<i a; 2_j / 2 superconducting state compared to the normal state is consistent with the 
experimental observations both at the zone edge and along the diagonal. However, the 
predicted peak widths in the superconducting state are narrower than those in the normal 
state, a narrowing which has not been observed experimentally. 

II. IMPURITY SCATTERING EFFECTS 

In this section we examine the effect of impurities on the spin susceptibility of a BCS 
d x 2_ y 2 state. We begin by including the self-energy effects and then consider the impurity 
vertex corrections. 

A. Self-energy corrections 

Neglecting the vertex corrections, the spin susceptibility in the superconducting state is 
given by the particle-hole bubble diagrams shown in Fig. 0. The bare vertex representing 
magnetic neutron scattering is understood to connect only electrons of opposite spin. Evalu- 
ating this diagram, continuing to real frequencies, and taking the imaginary part, we obtain 
the following formula for x"(q, tu): 

X "(q, W ) = -7T J j^- dJ [f(u + W ') - /(a/)] t^v(P + Q, " + «J)A N (p, u') 

+A A (p + q,u + u')A A (p,u;')}. (1) 
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Aw(p, w) and Aa(p, u;) are the spectral weights of the normal and anomalous electron prop- 
agators, respectively: 

A N (p,cu) = ImG(p,uj), (2a) 

7T 

A a { p ,uj) = — ImFfau), (2b) 

TV 

where 

G(p,u)= _. 9 Q l ip ~ (3a) 

and 

W 2 - £2 - A p 

The tilde symbol indicates inclusion of the impurity scattering self-energy corrections: 

u>— uj — S , (4a) 

e p = e p + S 3 , (4b) 

A P =A P + Ei. (4c) 

Here we are using e p = —2t(cosp x + cosp y ) — \i and A p = oA(T) (cosp x — cosp^). The 
parameter a is chosen such that the maximum value of the gap on the Fermi surface is A (T) . 
The effect of impurity scattering may now be included by allowing the electron self 
energy to include multiple scattering from the potentials V. Figure |3| shows these scattering 
processes for the dilute or non-crossing limit. Electrons are allowed to scatter multiply from 
an impurity to allow for arbitrarily strong scattering strengths. In the dilute limit, however, 
the electrons are assumed to interact with only one impurity at a time so the impurity 
interaction lines may not cross. The self energy within this approximation is given bya 

E (w) = T n T (uj), (5a) 

E 3 (w) = T n T 3 (lu), (5b) 
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FIG. 3. Diagrammatic representation of the calculation of the electron self energy due to im- 
purity scattering in the non-crossing or dilute limit. 



where 



T ^ = ?#&■ (fa » 



and 



T3M = cS-GgM " (6b) 

Here T N = ni/[irN(0)}, c = cot S , and 

G ° (W) = ^W) l (270^2 _^_ A 2' (7) 

where r^ is the impurity concentration, iV(0) is the normal phase density of states, and 8q is 
the scattering phase shift.ty The self-energy correction to the gap function Ei vanishes for 
a d-wave gap, and in the unitary limit, c = 0, only the E contribution remains. 

We can get an idea of the behavior of the electron self energy due to impurity scattering 
by looking at the quasiparticle decay rate 1/T imp {uj) implied by the self energy: 

l/r imp (uj) = -2ImE (w). (8) 



Figure £| shows 1/T imp (u) at T = 0.1T C for n« = 1%, n« = 0.3%, and n« = 0.1% in the 
resonant or unitary limit, c = 0. Here X (co>) nas been evaluated for a filling slightly below 
half filling < n >= 0.85 and a low temperature gap magnitude given by A(0) = 3T C , where 
T c = 0.05t, parameters which are also used in the calculations that follow. 

The unitary limit is of some interest in these systems as it has been invoked to explain the 
T 2 temperature dependence of the magnetic penetration depth seen in the thin film cuprate 
superconductors.^ In general, unitary limit scattering leads to the largest low frequency 
quasiparticle decay rate for a given impurity concentration. For very large quasiparticle 
decay rates, l/r imp (uj) ~ A(T), we might expect the impurity scattering will lead to a 
suppression of the superconducting gap itself, possibly leading to a reentrant behavior of 
the superconducting state. In the calculations which follow we use n« = 1% with c = 0, a 
choice of parameters such that 1/T imp (uj) is always smaller than A(0). 

Figure || shows the results of a calculation of x"(q, w ) with and without the inclusion of 
impurity self-energy effects. Each of the following figures shows x"(q, cj) evaluated at low 
frequency, u = 0.4T C , for q values which run along the solid line in Fig. [I], a line which passes 
through the normal state peaks in x"^ uS). The x-axis shows the distance Q in momentum 
space measured from the point Q 7 = {it — 5/2, tc — 5/2) where the zone diagonal (dashed 
line of Fig. [[]) intersects the solid line in Fig. [I], that is, q = (n — 5/2 + Q, n — 5/2 — Q). For 
the above filling, the normal state peaks are displaced from (tc, tc) by an amount 5 = 0.9ir. 
This value for 5 depends upon the details of the band structure. For a three-band Hubbard 
model with parameters adjusted to fit the La2Cu04 band structure it can be about twice 
as largeffl in agreement with experiments on La2_ x Sr x Cu04. Here for simplicity, as in some 
previous studies,0 we use a one-band near- neighbor form for e p . 

In Fig. |^ curves are shown comparing the normal, T = T c , and superconducting, T = 
0.1T C , states as well as the pure and dirty cases. For the normal state calculation, we would 
expect that the electron self energy should include contributions due to inelastic scattering 
by spin fluctuations (this is discussed further in Sec. fT|). Here instead, the self energy in 
the normal state has been modeled as a frequency and momentum independent constant, 
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FIG. 4. Quasiparticle relaxation rate l/rj mp due to impurity scattering in the unitary limit, 
c = 0. Results are shown for n,i = 1% (solid line), rii = 0.3% (dashed line) and n, = 0.1% (dotted 
line). 
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FIG. 5. x"(q, w) with and without the inclusion of impurity self-energy effects as a function of 
momentum. The momentum scan extends along the solid line in Fig. Ill perpendicular to the zone 
diagonal. Q = indicates the intersection with the zone diagonal. The vertical axis is in units of 
t . The curves represent the following: (solid) normal state, pure; (short dashes) normal state, 
dirty; (long dashes) superconducting state, pure; (dotted) superconducting state, dirty. 
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So = — ij, where 7 = T c /2. It is apparent from the figure that this self energy has only a 
modest effect on the shape of the x"(q, a>) peaks. Thus the use of such a simple model for 
the self energy should not interfere with interpretation of the experimental results. 

In the superconducting state the inelastic scattering due to spin fluctuations is suppressed 
by the opening of a superconducting gap. Thus impurity scattering effects will dominate at 
low temperatures and frequencies. Accordingly, the dirty superconducting state calculation 
incorporates an impurity scattering self energy calculated as shown in Eq. fl5|). Figure || 
shows the results for an impurity concentration ni mp = 1% with a scattering strength in the 
unitary limit, c = 0. 

As may be seen in Fig. |5], for the pure case calculation the zone edge peaks which occur in 
the normal state are greatly suppressed in the superconducting state and a peak in x"(q, w ) 
occurs instead along the zone diagonal. Inclusion of impurity scattering self-energy effects, 
however, works to restore the zone edge peaks while the zone diagonal peak is washed out. 
A lack of thermally available scattering states leads to the suppression of the zone edge 
peaks for the clean superconductor. In the dirty system, impurity scattering broadens the 
quasiparticle resonances in the electron propagators allowing access to otherwise thermally 
restricted scattering states and acting to restore the zone edge peaks which arise due to 
scattering between such states. 

B. Vertex corrections 

The above electron self energy includes all possible non-crossing impurity scattering con- 
tributions to the electron self energy. In order to include all possible non-crossing impurity 
scattering contributions to the spin susceptibility, we must also include vertex corrections 
due to impurity scattering. That is, we must include all processes where the electron and 
hole are allowed to scatter multiply from the same impurity. The form of these vertex cor- 
rections is shown in Fig. |^. The vertex corrected susceptibility bubble %(q, iuj m ) is given by 
a particle-hole bubble M(q, iu> m ,iuj n ) with a simple vertex at one end and a dressed vertex 
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FIG. 6. Diagrammatic representation of the vertex corrections to the spin susceptibility in the 
normal state due to impurity scattering in the non-crossing or dilute limit. 
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A(q, iu m ,iuj n ) at the other end. In the normal state this may be written as 

x(q,iuJm) = TJ2 M(q,iu> m ,iu) n ) A(q,iuj m ,iuj n ), (9) 

n 

where 

/d v 
——^G(p + q,iu m + iu n )G(p,iu n ). (10) 



(2k) 

The dressed vertex is given by the sum of a series of ladder diagrams where impurity scat- 
tering lines [represented by Y{iuj m ,ioj n )} are allowed to connect the particle and hole lines: 

A(q, iuj m , iu n ) = 1 + T(iu m , iu n ) M(q, iu m , iu n ) A(q, iu m , iu n ) . (11) 

Equations (jfty and (|TT|) are easily rearranged to yield 

Wq,»w m = r Li 7T- : — TUT : : — 7- ( 12 

Since the impurity scattering lines include any number of multiple scattering events, they 
are given simply by a product of two of the t-matrices used in the self energy calculation: 

r(iu m ,iu n ) = — % T (iuj m + iu n ) T (iu n ). (13) 

Since we are considering only the unitary scattering limit, i.e. c = 0, the T3 component of 
the scattering t-matrix vanishes. Thus T(iu m ,iu n ) contains contributions from T only. 

A complication which arises in the superconducting state is that the dressed vertex 
cannot be represented as a simple scalar function. Since the electron and hole propagators 
must be represented by both normal and anomalous propagators to account for electrons 
scattering into and out of the superconducting condensate, the dressed vertex is represented 
by a four component vector which accounts for all possible combinations of electron and 
hole lines going into and out of the vertex, as shown in Fig. [7|. The equations for the vertex 
corrected spin susceptibility 

X(q, iu m ) =T^2 M(q, iw m , iw n ) A(q, iu m , iu n ) (14) 
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FIG. 7. Components of the dressed vertex in the superconducting state. 
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and for the dressed vertex 

A(q, iu m , iu n ) = T + T(iu m , iu n ) M(q, iu m , iu n ) A(q, iu m , iu n ) (15) 

then become matrix equations where the vertices are 4-component vectors (overscored sym- 
bols) and the particle- hole bubbles are represented by 4x4 matrices (underscored symbols). 
Equation ( |15D represents four separate equations, one for each component of A(q, iu m , iu n ). 
The diagrammatic representation of the first of these is shown in Fig. [8]. The diagrammatic 
representations for the other three components are similar. The diagrams corresponding to 
the components of the matrix M_(q^,iu m ,iuj n ) are shown schematically in Fig. |9|. Explicitly, 
the components of M(q. iujm, iu)„). which is symmetric, are 

M n (q, iu m , iid n )= -M 22 (q, iu) m + 2iu) n , -iu) n ) = -M3 3 (q, -iu m - 2itn n , itn n ) 

/d?v 
——2 G(p + q, ito m + iu n ) G(p, iu n ), (16a) 

M 12 (q, iu m , iuo n )= M 13 (q, -icv m , iuj m + iw n ) = -M 24 (q, iuj m + 2iuj n , -iuo m - iuj n ) 

/d 2 p 
——2 G(p + q, iuJ m + w n ) F(p, w n ), (16b) 

/d?p 
——2 F(p + q, iu m + iu) n ) F(p, iu n ). (16c) 



Just as for the normal state, Eqs. (|I4"D and (|1"5|) may be rearranged, yielding 



X(q.,iuj m ) =TJ2 M(q, «*>m, «*>n) [l-T(iuj m ,iuj n ) M(q, zw m , zu; n )] *. (17) 

n 

As mentioned previously, the vertex representing magnetic neutron scattering connects 
only electrons of opposite spin. Thus to recover the quantity corresponding to the spin 
susceptibility measured by magnetic neutron scattering we must now add the components 
of x(q, iwm) which correspond to diagrams which begin and end with opposite spin electron 
lines: 

X(q,iu m ) = x n (q,^ m ) + x 41 (q,^ m )- (18) 

15 
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FIG. 8. Diagrammatic representation of the first component of the matrix equation, Eq. (|15D, 
for the dressed vertex in the superconducting state. 
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To avoid double counting, X 14 (q;2Co> m ) an d ^(1%) are not included. 

Figure |10| shows x"(q, w ) calculated with and without impurity vertex corrections. Here 
x(q, iuj m ) has been calculated for imaginary frequencies. The real frequency results were 
then obtained by analytic continuation using Pade approximants.liil In the normal state the 
vertex corrections do not have a dramatic effect on the zone edge peak structures. There 
is a noticeable decrease in x"(q, to) along the zone diagonal. In the superconducting state 
the impurity vertex corrections lead to a further enhancement of the peaks at the zone edge 
positions. These peaks are now comparable in size to those in the normal state. 

III. COULOMB INTERACTION EFFECTS 

In addition to impurity scattering, it is important to take into account the effects of the 
spin correlations produced by the Coulomb interaction. These have been shown to play an 
essential role in calculations of both the normal and the superconducting NMR responses.0 
Furthermore, as noted in the introduction, Coulomb spin susceptibility enhancement effects 
have been shown to maintain the neutron scattering peaks at their normal state positions 
along the Brillouin zone edge.&0 Just as in the NMR and neutron scattering studies, we 
treat the Coulomb interaction within an RPA approximation in which 

/ s Xo(q,w) . , 

x(q,w) = - — jt—, — v ( 19 ) 

Here, however, xo(q;^) is the spin susceptibility in the presence of impurities discussed in 
Sec.[n[ The interaction U represents a renormalized Coulomb interaction.113 Figure |ll| shows 
x"(q, oj) calculated for U = It. The Coulomb interaction enhances the normal state peak 
height by close to a factor of 60. This results in peak heights for the model which are of 
the same order of magnitude as the spin susceptibility peak heights reported by Mason and 
co- worker s .EltJfca 

As was mentioned in Sec. |TJ, the Coulomb interaction also gives rise to dynamic life- 
time effects through spin-fluctuation scattering. Previous calculations^ have shown that 
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FIG. 10. x"(q, u) with and without the inclusion of impurity vertex corrections as a function 
of momentum. The curves represent the following: (solid) normal state, vertex corrected; (short 
dashes) normal state, no vertex corrections; (long dashes) superconducting state, vertex corrected; 
(dotted) superconducting state, no vertex corrections. 
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FIG. 11. x"(q, u;) with Coulomb enhancement effects included in the normal (solid line) and 
superconducting states (dashed line) as a function of momentum. 
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scattering of electrons by spin fluctuations leads to the right order of magnitude for the 
quasiparticle relaxation rate in the normal state of the cuprates, i.e., t~ 1 (T c ) ~ T c . In 
the superconducting state, however, the gap in the electron excitation spectrum leads to 
a suppression of the low energy part of the spin fluctuation spectral weight and therefore 
to a suppression of the spin fluctuation scattering contribution to the quasiparticle relax- 
ation rate.£3 Previous calculationsE-3 for a model system similar to the one used here have 
shown that for T = 0.1T C and quasiparticle frequencies of order the temperature the quasi- 
particle relaxation rate due to spin fluctuation scattering is suppressed by several orders of 
magnitude compared to the value at T c . For higher frequencies, such as those used in the 
experiments of Mason and co-workers, this relaxation rate is still small compared to the 
impurity scattering rate for the concentration of impurities we have considered here.E-3 Thus 
the impurity scattering represents the dominant lifetime effect in the superconducting state 
results discussed here. In the normal state, on the other hand, the lifetime is determined by 
inelastic spin-fluctuation scattering. However, in the normal state, lifetime effects lead only 
to small changes in peak structure of x"(Q> to), as shown in Fig. [|. 

In order to directly compare the calculation including Coulomb enhancement effects to 
the experimental results we must average x"(q, tu) over a region of momentum space com- 
parable to the experimental resolution. Figure ^ shows x"(q, tu) averaged over a region in 
momentum space with dimensions 0.047T by 0.12-7T, the larger dimension being perpendicular 
to the momentum scan direction. These results show some qualitative agreement with the 
experimental data. At the zone edge and zone diagonal the magnitude of x"(q, w ) is reduced 
by a factor of about one half in the superconducting state compared to the normal state. 
This reduction is, however, not quite isotropic. The zone edge peak structures are noticeably 
narrower in the superconducting state compared to the normal state. 



20 



0.4 



0.3 - 



3 

cr 

E 



0.2 - 



0.1 - 








-0.2 -0.1 





Q/tt 



0.1 0.2 



FIG. 12. Momentum averaged x"(q, to) in the normal (solid line) and superconducting states 
(dashed line) as a function of momentum. 
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IV. CONCLUSIONS 

We have calculated the spin susceptibility at and below T c for a model d a ,2_ y 2-wave super- 
conductor with resonant impurity scattering and Coulomb correlations. Impurity scattering 
effects, including vertex corrections, act to restore the zone edge peaks in the spin suscep- 
tibility which are otherwise suppressed by a rf-wave superconducting gap. The predicted 
amount of suppression of the spin susceptibility in the superconducting state compared to 
the normal state is in qualitative agreement with the experimental results in that the pre- 
dicted neutron scattering intensity at the zone edge and along the diagonal is suppressed 
by about one half in the superconducting state compared to the normal state. The peak 
structure predicted by this model differs from the experimental results in that the peaks 
narrow appreciably at low temperatures compared to the normal state. 
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